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Abstract 

For the generalized statistical mechanics based on the Tsallis entropy, a variational perturbation 
approximation method with the principle of minimal sensitivity is developed by calculating the 
generalized free energy up to the third order in variational perturbation expansion. The approxi- 
mation up to the first order amounts to a variational approach which covers the variational method 
developed in Phys. Rev. Lett. 80, 218 (1998) by Lenzi et al, and the approximations up to higher 
orders can systematically improve variational result. As an illustrated example, the generalized 
free energy for a classical harmonic oscillator (considered in the Lenzi's joint work) are calculated 
up to the third order, and the resultant approximations up to the first, second, and third orders 
are numerically compared with the exact result. 
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I. INTRODUCTION 



Tsallis non-extensive statistical physics (TNESP) is one of theories for generahzing the 



fl. 



Its formahsm is based 



Boltzmann-Gibbs statistical mechanics and thermodynamics 
on the Tsallis entropy with a parameter g, the index of non-extensiveness for a system. 
It takes the conventional statistical mechanics as its special case oi q ^ 1. A distinct 
feature of it is the power-law distributions (g 7^ 1) instead of the exponential law in the 
conventional statistical mechanics. Since Tsallis' pioneered work in 1988 f^, the TNESP 
has greatly developed, and successfully been used for investigating many systems with a 
long-range interaction, long-time memory, or fractal structured space-time (l|, in which the 
non-extensive effect exists |l|. 

For a new physical theory, to develop basic and analytic approximation tools for the 
TNESP should be a fundamental task for developing and employing it. To finish this task is 
particularly necessary for the TNESP because to perform calculation with the Tsallis statis- 
tics is usually difficult due to the presence of the power-law distributions. Since the TNESP 
with g = 1 is reduced to the conventional statistical mechanics, it will be natural to gener- 
alize variational method and perturbational theory, the two basic approximation methods 
in the conventional statistical mechanics [sl], to the TNESP. In fact, some progress has been 
made in this direction. Early in 1993, employing concavity properties of Tsallis entropy, 
Plastino and Tsallis established a generalized Bogoliubov inequality and accordingly devel- 
oped a scheme of variational approximation for the free energy in the generalized canonical 
ensemble j^. Later in 1998, Lenzi et al. generalized the perturbation method by deriving 
approximate expressions of the generalized free energy (GFE) up to the second order and 
obtained, through analyzing the approximate expression up to the second order, another 
generalized Bogoliubov inequality which affords a different variational method J5[ . Recently, 



brementioned variational methods ICl , and 



y, 



Mendes et al. gave a comparative study on the a 
illustrated that the variational method in Ref. ^ gives a better approximation for g < 1 
and is easier to perform than the variational method in Ref. jj|, albeit the latter can lead 
to a better approximation than the former for the case of 1 < g < 2. 

Nevertheless, the usefulness of both variational and perturbation methods is limited. 

n 

The variational approximate scheme in Ref. |^ is valid only for the case of g < 2, and the 







variational approximate scheme in Ref. |^ works only for the case of g > 0. As for the 



perturbation method, the nth-order perturbation approximation can be used only when q > 
(1 — -). Moreover, as is well known in the conventional statistical mechanics, to control the 
approximate accuracy of variational methods is not straightforward, and the perturbation 
method is valid only for an exactly soluble system with a really small perturbation. Hence, 
for systematically improving variational method and extending the valid range of q, it is 
worthwhile developing new approximate approaches. 

In the conventional statistical mechanics and some other branches of physics, a variational 
perturbation idea [7[ which collects merits and overcomes drawbacks of both the perturba- 
tion and the variational methods has developed and now becomes a powerful tool called 
variational perturbation theory [8] (sometimes nominated other names.). The author or 



with his col 



^(Refs. labl 



aborators also developed some schemes for the variational perturbation theory 



have briefly introduced references on the schemes for variational perturbation 
theory). In the present paper, we intend to generalize the variational perturbation theory 
to the TNESP. In the generalized canonical ensemble, we will consider an expansion for the 
GFE analogous to that in Ref. p but with artificially introduced some auxiliary parameter 
/i and an expansion index e, and derive the expressions of the first four terms in the expan- 
sion series which will lead to approximations for the GFE up to the third order. From the 
truncated expressions of the GFE at the first, second and third orders (with taking e = 1), 
the auxiliary parameter n will be determined according to a principle of minimal sensitivity 
(PMS) jlO] at the first, second and third orders, respectively. The value of ft determined at 
some order entering the truncated expression of the GFE at the same order produces the 
approximation for the GFE up to the same order. One will see that the adoption of the 
PMS is crucial to the variational perturbation approximation (VPA) method. It is the PMS 
that guarantees the meaningfulness, non-perturbational nature and effectiveness of the VPA 
results. To illustrate the method, we will consider a classical harmonic oscillator, which was 
used in Ref. p| , and calculate the GFE up to the third order in the variational perturbation 
expansion. Numerical comparisons between the first-, second- and third-order and exact 
results are made and indicate that the VPA method provides a better approximation than 



the variational method in Ref. 

Nowadays, owing to the so-called "normalization" problem of ^. 
mostatistics, there have existed four versions for the TNESP j^, |ll 



with the usual constraint on the internal energy in Ref. [2] was used, for some special sys- 



a 



sallis non-extensive ther- 



12 



13l |. The first version 



terns, only a couple of times in the past, and the second version with the non-normalized 
constraint on the internal energy in Ref. [ll| has been intensively studied and used, and fur- 
thermore the aforementioned variational and perturbation methods were developed in this 
version. The third version with the normalized constraint on the internal energy in Ref. 12 1 
is a satisfactory version, but is very complicated for performing because there exist implied 
relations between relevant quantities, and so the fourth version with the non-normalized 
constraint on centered operators appears to unentangle the problem on the implied relations 
in the third version |l3|. Happily, the four versions can be easily derived from just any one 



of them 



12 



IJ]. Moreover, we notice that the second version has been applied to many 



systems with providing satisfactory theoretical and/or experimental results [l|. 



]J. Hence, 



llj l to perform the VPA method. 



here we will use the second version worked out in Ref. 

Next, the Tsallis statistical mechanics with the canonical ensemble and the perturbation 
and variational methods in Ref. [5| will briefly be introduced for the convenience of our later 
investigations. In Sect. HI, the VPA method will be stated and truncated expressions of the 
expansion series of the GFE in the second version in Ref. |ll| for a system will be derived up 
to the third order. Sect. IV will contribute to investigating the classical harmonic oscillator, 
and conclusions will be made in Sect. V. 

II. TSALLIS STATISTICAL MECHANICS 

For a system with a non-extensive index g, the Tsallis entropy is defined as 

S, = k ^_^=f\ (1) 

where k is the Boltzmann constant and W is the total number of microscopic possibilities pi 
of the system. Note that g is a real number and characterizes the degree of non-extensiveness. 
When g = 1, Eq.(l) leads to the usual entropy 5*1 = —kY^^^^pihipi. The three cases of 
g < 1, g = 1 and g > 1 characterize superextensivness, extensiveness and subextensiveness 
of systems, respectively. Based on the Tsallis entropy, a generalized equilibrium thermo- 
statistics can be established by making Sq extremal with appropriate constraints present. 
For the canonical ensemble, four possible candidates exist for the constraints and hence lead 
to the aforementioned four versions for Tsallis statistical mechanics. This section will use 



the following non- normalized constraint on the internal energy Ug |lll |: 

w 
Y.plE, = Ug, (2) 

where Ei is the ith eigenvalue of the Hamiltonian of the system. (The corresponding con- 
straint used in the first version of Tsallis statistical mechanics is X]j=i PiEi = Ug, and the one 
in the third version is the normalized one X]i=iP?-^«/ Si=iP? ~ ^q-) Thus, making Sq ex- 
tremal with present the constraint Eq.(2) and the normalization property of the probability 



yields the following power-law distribution 



ul 



Pi = PiiEi) = PiiP ) = zUF) ^^ 

with the generalized partition function 

w 
Z,(/?*) = ^[l-(l-g)/?*E,]V(i-^), (4) 

where (3* is the Lagrange multiplier associated with the internal-energy constraint, Eq.(2). 
In the case of g < 1, the summation in the distribution Eq.(2) will be cut off for those 
energy eigenvalues higher enough to give negative probabilities. It is easy to verify that 
when q —>■ 1, the power-law distribution Eq.(2) tends to the conventional exponential-law 
distribution. 

From Eqs.(2), (3) and (4), one can easily read 

,, ^ d Z]-^-^ ^ dlUgjZg) 

" dl3* l-q d(3* ' ^^ 

where the g-logarithm function lnq(a;) = {x^~'^ — 1)/(1 — q) is a generalization of the usual 
logarithm function. Through Legendre transform on lng{Zq), which depends on /?*, one can 
find the relation Sq = k{lng{Zg) + f3*Uq). Note that the relation between the generalized 
internal energy and Tsallis entropy takes the same form as in the conventional statistical 
mechanics. In fact, introducing t = l/{kf3*), and defining generalized thermodynamic func- 
tion in the same way as that in the conventional thermodynamics, one can find the same 
Legendre structure as in the conventional thermodynamics. For example, one can have the 
GFE 

Fg^Uq- tSg = Fq{n = -J.^^ = -^ 1^'^(^'^) • (6) 



Using the basic limit formula \iTa;j.^Q{l + ax)^ = e"", one can verify that the above generalized 
thermodynamic function and partition function are reduced to the usual ones in the limit 
of g — > 1. 

From the above, one can see that although the usual thermodynamical Legendre structure 
remains valid in Tsallis statistical mechanics, it is evident that the original calculation 
techniques can not be directly borrowed into the generalized theory. One has to design 
calculation techniques for the generalized theory. As was stated in the introduction, a 
united presentation of the perturbation and variational methods for the TNESP has been 
given in Ref. p]. Next, for convenience of later contrast, we give a brief introduction on 
them in the notations here. 

Assume that the Hamiltonian can be written as 

H = Ho + XHj, (7) 

where, Hq is the Hamiltonian of a soluble model, XHj is small enough so that it can be 
considered as a perturbation on Hq, and A is the perturbation parameter. The GFE of the 
system with H, Fg{X), is a function of A and can be expanded as 

i^.(A) = Fr + AF« + ^Ff + ... (8) 

with Fq = g| ' |a=o, Fq = g^2 |a=0; " " " • Eq.(8) is a perturbation expansion series of 
Fq{X) in powers of A, and to truncate it at some order, for example nth order, can provide 
an approximation to Fq{X) up to the nth order. In Ref. IS], approximate expressions for 
Fq{X) was derived up to the second order. Through analyzing the approximation expression 
of Fq{X) up to the second order, Ref. p obtained the generalized Bogoliubov inequality 

w 

v-^ m^ (0) - (0) 

(0) 

where the symbol O is the operator corresponding to the physical quantity O and \i > (and 

(0) 

< i\) is the ith eigenstate for Hq. This inequality stems from the approximation expression 
of Fq{X) up to the first order in the perturbation expansion, and affords a variational method. 
According to Ref. p|, the approximate expression of Fq{X) up to the nth order is valid 
only for the case of g > 1 — - because the derivation of the approximate expressions requires 
the interchangeability between the sum over the quantum numbers and the derivative with 



respect to A. As was stated below Eq.(4), when q < 1, the sum in Eq.(4) and, accordingly, 
Eq.(6) has a cutoff and so its upper limit depends upon A, leading to the above limitation 
on the resultant expressions in Ref. |5|. Therefore, the variational method in Ref. J5| can be 
used only for g > 0. 

We stop here for introducing the perturbation and variational methods in Ref. [5|. Next, 
we turn to our investigation in the present paper. 

III. VARIATIONAL PERTURBATION APPROXIMATION SCHEME 

For simplicity of the notation, k = 1 will be taken from now on. 

To develop the VPA method, we consider calculating the GFE Fg (Of course, the method 
can be used to calculate other generalized thermodynamic potentials). We begin with mod- 
ifying the Hamiltonian H. Firstly, differently from Eq.(7) in the perturbation method, we 
add the zero term Hq{^) — Hq{^) to H and write it as 

H = Hoifi) + Hj{fi) (10) 

with Hj = H — Holfi) and /i an auxiliary parameter. In Eq.(lO), Hq{^) is the Hamiltonian 
of some exactly soluble model which is originally not included in if, and the exact solubility 
of the model is not affected by values of the auxiliary parameter /x. That is to say, for a 
value of /i, the energy eigenvalues En' (non-degenerate, for simplicity), eigenstates \n >, 
the probability Pq{En ), the generalized partition function Zg^ and the GFE F^ for the 
model with Hq{^) are exactly known. Secondly, an artificial parameter e is inserted as a 
factor before Hj{fi) in Eq.(lO), that is, H in Eq.(lO) is now modified as 

H, = Ho{ij.) + eHi{^). (11) 

Thus, we have a new system with the Hamiltonian H^ which becomes Ho{fi) and H in the 
cases of e = and e = 1, respectively. Note that different from H in Eq.(lO), H^ depends 
on fi for € ^ 1. The above modification provides a possibility that we calculate Fg for the 
original system through considering the new system. 

For the system with H^, suppose that the set of energy eigenvalues are the spectra En, e, 
and the corresponding eigenf unctions is \n,e >. When the system is in thermal contact with 
a reservoir, i.e., when we adopt a generalized canonical ensemble, we have the probability 



distributions pn{En,e), and can calculate the GFE Fq{e) in the presence of the constraint 
Eq.(2) according to the definition and expression Eq.(6). In analogy to Eq.(8), one can 
expand Fg{e) as the following series 

F,ie) = Ff + 6F(^) + ^Ff + ^Ff + ■ ■ ■ (12) 

with Fq = g^f^ I _ and Fq = Fq(0). This is only a formal expansion, and the parameter 
e acts as an expansion order index and needn't be small. Taking e = 1 in the last equation, 
we obtain an expansion series of Fq for the system with H which is independent of n if the 
resultant series is not truncated. 

Now we have a simple discussion on the expansion Eq.(12). Of course, if Hq{^) is naively 
one part of the original H and Hj can be regarded as a perturbation on Hq{^) (in this 
case, /i is one of system parameters in H), then Eq.(12) with e = 1 is just the perturbation 
expansion series of F^, Eq.(8). If ifo(/^) is not a naive part of if , as we assumed in Eq.(lO), the 
truncated expression at the first order in e can lead to a variational method, and substituting 
the value /x from the variational method into Eq.(12) with e = 1 will produce an expansion 
series of Fq around the variational result. In this case, if the resultant variational method can 
provide a good approximation result for F^, then it can be presumed that the expansion of Fq 
around the variational result maybe afford a good non-perturbational approximation which 
will improve variational results. This is the original version of the variational perturbation 
theory, but the expansion series is possibly divergent [8] . If the variational method cannot 
provide a good approximation result for Fg, the expansion of Fq around the variational result 
will make no senses. That is to say, in general, the above expansion series Eq.(12) seems to 
be useless. 

However, as was mentioned in Sect. I, once the PMS is employed to single out an appro- 
priate value of /i, Eq.(12) can afford a presumably convergent, systematic non-perturbation 

n 

approximation method, a modern version of the variational perturbation theory 8]. In 
Eq.(12), truncating the expansion series at some order in e and then letting e = 1, the 
resultant truncated expression is a function of the parameter /x, albeit Eq.(12) with e = 1 is 
independent of /x. In principle, the value of the truncated expression can possibly take any 
value with variance of /i, including the exact value of Fq. Since Fq is a constant in the space 
of n (the exact Fq is independent of /i), it should be believable that a requirement of the 
truncated result varying most slowly with the parameter /x can likely make the truncated 

8 



result provide a most reliable approximate result for the exact Fq. Employing this require- 
ment to determine the parameter /x is the main spirit of the so called PMS 10]. Generally, 
the curve of a function varies more slowly near its extremum than in its slope part, and so 
a simple realization of the PMS is to require the first derivative of the truncated expression 
with respect to /i to be zero. Thus, substituting the appropriate value of /i determined at the 
truncated order into the corresponding truncated expression of Fq will provide a reasonable 
and reliable approximation for Fq up to the same order. Such a truncated result at some 
order is the approximate result of Fq up to the same order in the VPA method. 

In brief, the VPA method consists in only two crucial steps: one is to formally expand 
the quantity Fq in consideration with an adjustable parameter /x entered, and the other is to 
determine the value of /i from the truncated expression of Fq according to the PMS. Now, 
we continue to perform the first step, derive the first four terms in Eq.(12) (the derivations 
here are similar to those in Ref. J5|), and simultaneously indicate the VPA procedures up to 
the first, second and third orders, respectively. 

The first term in Eq.(12) is straightforwardly 

with 

w 

zun = E[i - (1 - <iWE^:^"^'-^^- (14) 

n=l 

For deriving the other three terms, eigenstates of H^ are needed, and so we first consider 
the eigenequation, Hf\n,e >= En.e\n,e >. By mimicing Rayleigh-Schrodinger perturbation 
theory, after substituting the expansions 

|n, e >= \n > +e\n > +e \n > -\ — ■ (15) 

and 

K,. = i?i°) + 6E«+6^E(^) + --- (16) 

into the eigenequation, one can formally have 



m^n 



F,^ ' - E 



m 



\n > 



2^ ( p(0) p(0)wp(0) 



-n/^ jYinJ^I, nn 



mj^n 



El 



{0)^ 



n(0) 






(0)^2 



m > 



(^, 



/, mnj 



2 Z^ /^(O) _ ^ 



\n > 



(18) 



(0) . (0) 

and so on. Here, Hjrnn =< m\Hj{fi)\n > . By the way, the expansion expressions, Eqs.(15) 
and (16), together with the PMS can lead to a VPA method to solve Schrodinger equation 

y(PRD)y. 

For the second term of Eq.(12), Fq = g^^' | can be calculated as 



ir(i) 

9 



d 1 - ^,^7 
9e/5*(l-g) 



e=0 



VeE[l-(l-^)^*^" 



.dE„ 



l-q . 



dt 



(19) 



e=0 



with Zq^^ the generalized partition function for the system with H^. From Hellmann- 
Feynman theorem —^ 



< e, n|^^|n, e > lid, we have 



g 



Zq^lY^^l - {I - q)f3*Er,,,]^^ < e,n\Hj\n,e > 



e=0 



= J2pl{E'jHi,nrr=<Hj>'^, (20) 

n 

which takes the same form as Eq.(5) in Ref. [5|. Thus, the truncated expression for Fq at 
the first order, F^{fi), is 



FUf^)=E' + FW 



yl-q 



(21) 



(3*{l-q) 

The right hand side of last equation has the same form as that in the right hand side of the 
inequahty, Eq.(9). F^{^) is a function of fi and a different value of /i produces a different 
approximation for Fg which maybe have nothing to do with the exact value of Fq. However, 
if the value of fi in Eq.(21) is chosen from roots of the following condition 



djj, 







(22) 



according to the PMS, then Fg{fi) in Eq.(21) would produce a most reliable approximation 
for Fq up to the first order. Since Fq has the same form as the right hand side of the 
inequality, Eq.(9), the VPA procedure up to the first order for the case of g > is consistent 
with the variational method based on the Bogoliubov inequality in Ref. [3], and so the 
appropriate root of Eq.(22) for n makes FUii) minimal. For the case of g < 0, an example 



10 



in next section shows that the appropriate root of Eq.(22) for /i also makes F^{fi) minimal. 
This suggests that the VPA procedure up to the first order is generally a variational method 
which provides a reliable upper limit for the exact Fg. 

In order to give the truncated expression for the GFE at the second order F^\fi) = 
F^(n) + Fq /2!, we calculate Fg as follows. A straightforward differentiation gives 



i?(2) 



d'FJe) 



de^ 



p*qz:>j 



e=Q 



;< Hj >y 



~Y.^P,{E^^)f''-\H,,^,y +Y.{Pg{E^^)) 



d < e,n\Hi\n, e > 



de 



(23) 



e=0 



From Eq.(15), we have 



de I 



(1) 



^= \n > . Consequently, Eq.(23) becomes 



Fj,'^ = P*QZIo 



[<Hi>lf-Y.{pM^)f''\Hj, 



+'^Y.^P,{E^^)YY. 



\H 



I, nm\ 



m^n 



p(0) _ p(0) ' 



(24) 



Formally, the right hand side of Eq.(24) can be rewritten as Eq.(6) in Ref. p|. According to 
the PMS, using the reasonable root of the condition 



dFl\fi) 







(25) 



as the value of /i, one can give the approximation result for Fg up to the second order from 



F (/i). The approximation result up to the second order would improve the approximation 



result up to the first order. 



.(3) 



Similarly, one can calculate Fq . In the calculation, we need the additional relation 
dJi'^ 1^-0= 2|ri >*^^\ which can be easily obtained from Eqs.(15), (17) and (18). Thus, Fg 
can be written as 



f(3) 
1 



P*\{q + l)Z, 



2q-2. 
q, > 



<Ht> 



0n3 



3/3*^g^<V'E(P^(^i°^))''"(^^."")' <Hi> 



+P*\{2q-l)Zl%'Y.ME^^)r-\H,, 



+6(3*qZ'gJ Yl 



H JnlH I JmH I ^nm 



I n |2 



(0) ZT.(0) 



vn^n 



E; 



E' 



+6 2^(Pg(^n )f l^l^ . (0) _ p(0)wp(0) _ p(0)^ A^ . p{0) 



\TT 12 IT 



l^r. 



r;(0)N2 



■ (26) 
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Then one can have the truncated expression for Fg at the third order , F^^^(/i) = F^\fi) + 
^Fq . Furthermore, using the reasonable root of the condition 

^^ = (27) 

as the value of fi, one can give the approximation result for Fg up to the third order from 

In the same way, one can consider VPA results for Fg up to higher orders. Obviously, 
the VPA results up to various orders constitute a sequence of the approximation results 
of Fg. Owing to the employment of the PMS, the values of /i determined depend upon 
the truncated orders, and the value of /x at a given order is different from values of /i at 
other orders. Presumably, it is the dependence of fi upon the truncated order that makes the 
sequence of the VPA results of Fg up to various orders converge to the exact Fg. Although we 
cannot verify this convergent property of the variational perturbation theory, investigations 
in other fields have provided such a few examples J8|, |l7[ . Additionally, in the example of 
next section, generally, the VPA result up to the third order approaches more closely the 
exact Fg than the approximations up to the lower orders. 

The second crucial step plays a vital role in the VPA method, and one have to perform it 
very carefully. Sometimes, Eqs.(22), (25) and (27) can not provide appropriate roots for /i, 
and in those cases one should render the second derivative of the truncated expressions for 
Fg with respect to fi zero to determine /x. This can be understandable because the curve of 
a function varies more slowly near its knee than in its slope part (As a matter of fact, when 
there exist both an extremum and a knee, one should analyze behaviors of the truncated 
expressions as a function of /i near the extremum and knee to determine fi according to the 
PMS). In case the second derivative condition can not provide an appropriate root for /i, too, 
then the VPA method cannot provide an approximation result for Fg up to the truncated 
order in consideration, and one should further consider VPA result up to the next order. 
Furthermore, when there exist multi-roots of the vanishing derivatives for fi, one should 
choose to adopt the root near which the truncated result varies most slowly with /i. Next 
section, we will show how to determine /x according to the PMS in a concrete example. 

In concluding this section, we intend to emphasize one point. In the above calculations, 
the interchange between the differentiation with respect to e and the sum over eigenstates was 
involved. As was mentioned in last section, the interchangeability between the differentiation 

12 



and the sum leads to a limitation for the variational and perturbation methods in Ref. J5|. 
Nevertheless, for the VPA method developed here, the relevant calculations are only formal 
calculations and the resultant expressions are independent of e, and presumablj^we suggest 
that the VPA method needn't suffer those limitations on the methods in Ref. p. At least, 
this point yields no problem in the example of next section. 

IV. A CLASSICAL HARMONIC OSCILLATOR 

For a one- dimensional classical harmonic oscillator with mass M and an angular frequency 
uj, the Hamiltonian is if = p^ / {2M) + M uj"^ x^ / 2 with x the coordinate and p the momentum. 
In a manner analogous to the one used for the quantum oscillator, we can associate with the 
oscillator a continuous energy spectra En = SoU, where 60 is an arbitrary positive constant 
with the dimension of energy and n any positive real number including jl^. Then, from 
Eq.(4), the generalized partition function can be written and calculated as 

Z, = l 11 -(l-,)3-„y. /(-,.,„=___ (28) 

with N ^ 00 for the case of g > 1 and A^ = 1/[(1 — q)f3*6o] for the case of g < 1 (In this 
section, we suppose that (3* > 0.). The right hand side of Eq.(28) is valid only for q < 2 and 
Zq is divergent for q > 2. So, Eq.(6) leads to the GFE for the classical harmonic oscillator 

^' ^ ~{l-q)(3*^^{2-q)P*6/'' ~ ^]- ^^^^ 

Last equation with 5o = huj is just Eq.(14) in Ref. ||5|. 

On the other hand, for a particle with mass M which moves in a one-dimensional box, 
the Hamiltonian Ho{L) = p^/{2M) + Vo with V"o = for |2;| < | and Vq ^ 00 for |a;| > f. 
In analogy to what we did for the oscillator, we can associate with the particle a continuous 
energy spectra En = 5ln^ / {2M L?') ^ where 5^ is an arbitrary positive constant with the 
dimension of action and n any positive real number including 0. So the generalized partition 
function for the particle in the box can be written and calculated as 



Z,o= [1 - (1 - ciW^]'/^'-'^dn = I ''^ V 0^!)^^^^' -^^' '^' (30) 

with A^ — > 00 for the case of g > 1 and A^ = La/2M/((1 — q)P*)/6b for the case of g < 1. In 
Eq.(30), B{x, y) is the Beta function. Note that Zq^ is convergent only for the case of g < 3. 

13 



In the calculation of the above equation, the formulae 8.380(1) and 8.380(3) in Ref. |l8| were 
employed. Owing to B{x, y) = T{x)T{y)/T{x + y) with r(a;) the gamma function (8.384(1) 
in Ref. yj), Eq.(30) with 5^ = hn is nothing but Eq.(ll) in Ref. p. From Eqs.(6) and (30), 
the GFE for the classical particle in the box is easily calculated as 



with 



c, 



2Mtt ^(l^) \l-q 



Sb(3-q)\ /3*(l-'?)r(^ + i)' 



g < 1 






Since the above two systems are exactly solved, Ref. J5[ employed them to illustrate the 
variational method there. In the present paper, we will use them to illustrate the VPA 
method. That is, by regarding the Hamiltonian of the classical particle in the box as Hq{^) 
in eq.(lO) and using the VPA scheme in last section, we will calculate the GFE for the 
classical harmonic oscillator up to the third order and then make comparisons among the 
results up to the various orders and the exact result for the oscillator. Next, the present 
section will be divided into three subsections. In subsection A, the truncated variational 
perturbation expressions for the classical harmonic oscillator will be given up to the third 
order . The values of /i in the truncated expressions will be determined according to the PMS 
in subsection B, and the approximated values of the GFE will be calculated and compared 
with the exact result in subsection C. 

A. Truncated Expressions 

Take the width L of the box as the adjustable parameter /i in our VPA scheme. Then, 
Hj{fi) in last section is now 

l' Mui'^x^ II ^ L 

l^^-oo, |a:|>f ^ ^ 

For the classical particle in the box, the speed Vn = n6h/{ML) corresponds to the particle 
energy En ■ In a state with the energy En , the particle moves in the box back and forth at 
the constant speed Vn, and so the probability of finding the particle in one of the directions 
near x (|x| < L/2) is dx/{vnT) with T = 2L/vn (except for n = 0), but zero out of the box. 
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Thus, in the present case, the matrix elements Hj^nm appeared in Eqs.(20), (24) and (26) 
can be straightforwardly calculated as 



H 



I,nm ^nm^ 



^/2 Miu^x"^ dx 



6„ 



Mu'L 



2t2 



L/2 



2 2L "" 24 



(33) 



where Snm = for m 7^ n and (5„„ = 1. 

Now we are at the position to calculate the truncated expressions for Fq of the classical 
harmonic oscillator up to the third order from Eqs.(13),(20),(24) and (26). The zeroth-order 
expression is Eq.(31). At the first order, one can have, from Eq.(20), 



n ^ 



2ML2- 



Mlj^L 



2 t2 



iZ,,o] 



24 

Mcj^L^ 3 - g 
24 2~ 



■1 1 



L / 2M p _ 

L I 2M D/'l 9+1 



g < 1 



a. 



(34) 



The right hand side of last equation is identical to Eq.(12) in Ref. 
At the second order in e, Eq.(24) leads to 



9 



^*qzij 



[<Hj> 



0n2 



J2iP.{E^:')r'-\H,,^:^) 



9,0 



/?*gZ, 



2, ,4r4 



M^cjU 



'^ ^ (24)2 (^' 



^9,0j 



N 



[l-il-lW^r-''^^'-''dn 



q < 1 



■\^2 4r4 ( — / 2Af p / 1 g N 

^.MW(g-3)(l-g)2^^^ (35) 



(24)^ 



4 
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Finally, we consider the third order. Eq.(26), in the present case, is reduced to 

n 



P*\{<i + l)ZT,\FJ^'f - ^P'\'Z;X^^^F^' j\^ - (1 - <i)P 



2X2 



n^Sz . zg-i 



2ML2 



+/3*^g(2g - l)Z^;^, 






[1 - (1 - qW 



r?bl . ag-2 



2ML^ 



] 1-9 (in 



r^g(g+lXV'(i^i'^)'-3/9 



^^'?,o (24)2 ^'^ \_L 



2M 



-5 



1 9 



(l-g)/3*^V2' 1-g^' 



g < 1 



2<56 V (g-i)/3* 



^^(iJfEly), g>i 



+/3*^g(2g - \)Z-% 



(24)= 



2M 



1 2g-l' 



9-0 I'O/I^S 1 r / 2M D/'l 5g-3 



256 V («-l)/3*""^2' 2{g-l)^' 



5( 



g > 1 



^,2 M3o;6L6(l + g)(3-g)(g-l)3 ^ 
^ (24)3 8 ^ 



(36) 



Thus, collecting Eqs.(31),(34),(35) and (36), one can get F/(L), F/^(L) and F^^^{L), the 
truncated expressions for Fg of the classical oscillator at the first, second and third orders. 
When q < 1, the truncated expression for Fg at the third order is 



q,sup 



(L) 



1 



L 



2M71 



r(T^: 



P*{1 - q) V^6(3 - g) V /9*(1 - g) r(^ + \] 



/?*(!- g)(3-g) 



24 



-P* 



(24)^ 



+ P 



*2 ' 



(24)^ 



(g + l)(g-l)= 



(37) 



and when 3 > g > 1, it is 



:pin 

q,sub 



(L) 



L_ n^£Mizll.i-.V" 



P*{l-q) \'26,V P*{q-1) r(-ij) 



/?*(!- g) 



^Ma^^3-g) _ ^, M^cu^L^ (3 - g)(l - g)^ 



+P" 



24 2 " (24)2 

^^M^u^L^ (g + l)(g - 1)^(3 - g) 



(24)^ 



(38) 



For the case of g > 3, F^jJ^^{L) in Eq.(38) is divergent. In Eqs.(37) and (38), all the terms in 
which the powers in uj"^ are lower than the second and the third powers give the truncated 
expressions for Fg at the first and second orders , respectively. Next subsection, we will 
determine L in them according to the PMS to get approximations for Fg. 
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In the following subsections, when a numerical calculation is performed, we will take 
M = l,uj = l,(5b= 1/2, and we will consider only the case of g < 2 which is the convergent 



range of Fq for the oscillator. [5^ = 1/2 corresponds to the case h = 1 in Ref 



.0.; 



B. Determining L 

L must be positively real. 

Up to the first order in e, for the cases of both superextensiveness and subextensiveness, 
Eq.(22) leads to the equation 



■^ + Mu\-3 + qrL' 



0. 



(39) 



To see whether Fq{L) is most insensitive to L near the positive root of Eq.(39) or not, one 
can analyze varying property of F^{L) as a function of L. For illustration, the function 
FHL) is depicted for (3* = 1 and some values of q in Fig.l. Fig.l indicates that the positive 




FIG. 1: The dependence of Fg{L) upon L for some values of q and at j3* = 1. The solid, dashed, 
dot-dashed, dotted and long-dashed curves from right to left correspond to q = 1.5,0.5, —0.5, —1.5 
and —2.5, respectively. Every curve in this figure has a minimum at a positive root for L. 

root of Eq.(39) for L makes FHL) reach a minimum and simultaneously be most insensitive 
to L. So, we should choose this positive root as the value of L 

4.3 



a/2 



3-q'Mu^p 



2R*' 



(40) 
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0. 



which is identical to Eq.(13) in Ref. J5|- The minimum F^{U) provides an upper hmit for 

Fq. Thus, the VPA up to the first order is really a variational method and it, in the case of 

o > 0, is consistent with the variational method in Ref. l5| which is based on the generalized 

fl 
Bogoliubov inequality, Eq.(9) in Ref. jq]. 

Up to the second order , for the cases of both superextensiveness and subextensiveness, 

Eq.(25) leads to the equation 



4608/?*^' - 96L'M{q - 3y/3*-^uj' + L'^M^q - 1)^(15 - 8g + g')w^ = . 



(41) 



Last equation produces four roots for L. The roots are real only when —1.96239 < q < 
1.62222 and 3 < g < 5.34017, and two of them are positive as well as the other roots are 
negative. To choose an appropriate root as L, F^\L), as a function of L, is representatively 
illustrated for (3* = 1 and some values of q in Fig. 2. In Fig. 2, horizontal lines represent exact 
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FIG. 2: The dependence of F^ (L) upon L for some values of q and at (3* = 1. The dot-dashed 
and dotted curves are for q = 1.25 and 0.1, respectively, and the solid and dashed curves in the 
image are for q = —0.5 and —1.5, respectively. Horizontal lines represent exact values of Fq, and 
every pair of curve and horizontal line with an identical line-type is for the same value of q. In the 
figure, every curve has a minimum and a maximum. 

values of Fq, and every pair of curve and horizontal line with an identical line type is for 
the same value of q. From curves in Fig. 2, one can see that the smaller positive root makes 
F^\L) minimal whereas the larger positive root makes F^\L) maximal. Hence, F^^{L) is 
more insensitive to L near the smaller positive root than near the other positive root, and 
so the smaller root should be chosen as the value of L, L^' , for —1.96239 < q < 1.62222. 
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For the other values of g, Eq.(25) cannot produce reasonable values for L. In this case, 



the vanishing requirement for the second derivative of F^ (L) with respect to L should be 
considered and yields the following equation for L 



q{AQ08p*-' - 20161^ Mp*-'iu^ + 167LW^cj^) = 



(42) 



which is valid for the cases of both superextensiveness and subextensiveness. Eq.(42) 
produces four roots for L. There exist real roots only when —3.02227 < g < 3 and 
3.95406 < q < 5.601, and for g < 0, two of them are positive as well as the other roots 
are negative, and for g > 0, only one positive root exists. To choose an appropriate root 
as L for the case of —3.02227 < g < —1.96239, F^^^L), as a function of L, is typically 
illustrated for f3* = 1 and g = —2 in Fig. 3. For completeness, we also depict F^^{L) for 
g = 1.8 (which is involved in the range of 1.62222 < g < 2) in Fig.3 as an image. In this 
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FIG. 3: Similar to Fig. 2, and the solid curve is for q = —2 and the dashed curves in the image is 
for q = 1.8. The solid curve has two knees and the dashed line in the image has only one knee. 

figure, the solid curve for g = —2 has two knees which correspond to the two positive roots 
of Eq.(42) for L: 1.13365 and 1.94787. Obviously, F^^{L) is more insensitive to L near the 
larger positive root than near the other positive root. The same situation occurs for other 
negative values of g (—3.02227 < g < —1.96239), and so the larger positive root of Eq.(42) 
should be chosen as the value of L, L^\ for —3.02227 < g < —1.96239. Furthermore, as is 
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illustrated in the image of Fig.3, for the case of 2 > g > 1.62222, F^ (L) is most insensitive 
to L near the positive root of Eq.(42), and so the positive root of Eq.(42) should be chosen 
as the value of L, L^^ , for 1.62222 < q < 2. By the way, for the case of g < —3.02227, there 
exists not any appropriate root of Eq.(41) or Eq.(42) for L and so the VPA method cannot 
provide an approximation for Fg up to the second order in this case. 

Finally, we determine L so as to approximate Fg up to the third order. For the cases of 
both superextensiveness and subextensiveness, Eq.(27) leads to the equation 

-663552/?*-^ + 13824L2M(-3 + q)^(3*~^uj^ - U4L^M^{-1 + g)2(15 - 8q + q'^)f3*-^uj^ 
+L^M%-1 + q)%21 + llq - 9q^ + q^)u^ = (43) 

Last equation produces six roots for L. Among the six roots, there exist three positive 
roots for —3 < q < —1, two positive roots for — 1 < g < 1 and only one positive root for 
q < —3, 1 < g < 3 and g > 7. In Fig.4, F^^^L), as a function of L, is typically depicted for 
some values of g and /?*. In Fig.4, the dotted curve is for {q,P*} = {—2.9, 1}, and has two 

irlll 




FIG. 4: The dependence of F^^^L) upon L for some values of q and (3*. The 
solid, short-dashed, dot-dashed, dotted and long-dashed curves correspond to {q', /?*} = 
{1.5, 10}, {0.5, 10}, {-0.5, 1.9608}, {-2.9,1} and {-3.5,1}, respectively. The solid and the long- 
dashed curves have a minimum, the short-dashed and the dot-dashed curves have a minimum and 
a maximum, and the dotted curve has two minima and one maximum. 



minima and a maximum, which corresponds to the case of having three positive roots. In 



this case, it is evident, from Fig.4, that F (L) is most insensitive to L near the smallest 
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positive root of Eq.(43) (the curvature of the curve at the smallest positive root is smaller 
than that at the largest positive root), and so the smallest positive root should be chosen 
as the value of L, L^^^ , for — 3 < g < —1. The short-dashed and the dot-dashed curves are 
for {q,l3*} = {0.5, 10} and {—0.5, 1.9608} and have a minimum at the smaller positive root 
and a maximum at the other positive root. So, for the case of having two positive roots, the 
smaller positive root is appropriate for L when — 1 < g < 1. As for the cases of g < —3 and 
1 < g < 2, the solid and the long-dashed curves indicate that F^^\L) reaches an extremum, 
the minimum at the positive root, and so it can be taken as the value of L, L^^^ . 

According to the PMS, we have determined L for all cases which we are interested in. In 
the above, analytical and numerical discussions are done with Mathematica, and we do not 
list all expressions of roots of Eqs.(41), (42) and (43) because they are too lengthy. Note that 
although L is determined as L^ ,L^^ and L^^' for approximating Fq up to the first, second 
and third orders, respectively, it is not meant that L is being expanded as a series. 

C. Generalized Free energy and Comparisons 

From analysis and results in last subsection, employing Mathematica package, one can get, 
for various ranges of g, the expressions oi F^{L = L^), F^^{L = U^) and F^^^{L = L^^^), the 
approximations of the GFE Fg up to the first, second and third orders. Regarding them as 
functions oit {t = 1//5*) and the non-extensiveness index g, respectively, we can numerically 
calculate and compare them with the exact Fg, Eq.(29). 

The exact free energy, Fg in Eq.(29), reaches a maximum at t = 5o{2 — qyl^i^^^ for any 
given value of g < 2. As a function of g, Fg in Eq.(29) has a maximum for any t < 0.17 or 
so, and, otherwise, reaches first a maximum and then a minimum when g increases up to 
g = 2. 

In Ref. pI, the approximate GFE from variational method there was considered for < 
g < 2, and when t increases or when g decreases, the discrepancies between the approximate 
values and the exact values become more and more evident. Here, we consider the VPA GFE 
up to the first order, F^{L = L^), for all values of g < 2, which with < g < 2 is identical 
to that in Ref. jq. Basically, Fg{L = L^), as a function of t or g, mimics the feature of the 
exact free energy, and for a given value of t, when g approaches 2 or is sufficiently negative, 
the error AF = Fg{L = L^) — Fg is very small whereas AF is not small for intermediate 
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values of q. Interestingly, the dependent feature of AF upon q is similar to that of curvature 
of F^{L) at L = L^ , F" = ^^2 — \l=l'^ upon q (For subextensiveness case, when t is small, 
there exists not such a similarity.). For an illustration of this similarity, in Fig. 5, we depict 
the dependence of AF (dashed curve) and F" (solid curve) upon g at t = 3. In Fig.5, the 

AFF" 

1200 




FIG. 5: The dependence of AF (dashed curve) and F" (solid curve) upon q at t = 3. 

image is drawn for the range of— 12<g<0, and evidently indicates that while the solid 
curve reaches the maximum, the dashed curve also gets to its maximum. This similarity 
can be understood from the PMS. Smaller the curvature ai L = L^ is, more slowly F^[L) 
varies near L = L' , and so, closer F^{L = L^) approaches Fg according to the spirit of the 
PMS, yielding the similarity. 

The approximation of Fq up to the second order, F^^{L = L^^), improves the approx- 
imation of Fg up to the first order, F^{L = L^), the variational result for the range of 
—3.02 < q < 2 where F^\L = L^^) makes senses, and the approximation of Fq up to 
the third order, F^^\L = L^^^), generally further improves the variational result. For a 
comparison and an illustration, the exact and various approximate results are shown at 
t = 1 in Fig. 6. In Fig.6, the solid, short-dashed, dot-dashed and dotted curves are Fq, 
F^{L = L^), F^^L = U^) and F^^^L = L^^^), respectively (note that the dot-dashed 
curve, F^\L = L^^), interrupts at g = —3.02 or so, and almost coincides with the dotted 
curve, F^^\L = L^^^), when g > — 1 or so). Additionally, the long-dashed curve in Fig.6 
is Fg^\L = L^), and obviously it is a bad approximation for Fq when q > —2 or so, albeit 
F^\L = L^) can produce as a good approximation for Fq as F^^[L = L^^) does. Fig.6 
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FIG. 6: The dependence of Fg and its approximations up to various orders upon q at t = 1. 
The sohd, short-dashed, dot-dashed and dotted curves are the exact Fg and its approximations 
Fg{L = L^), Fg\L = L^^) and Fg^\L = L^^^), respectively. The dot-dashed curve interrupts at 
q = —3.02 or so, and almost coincides with the dotted curve when q > —1 or so. The long-dashed 
curve is Fg^\L = L ), and goes up quickly when q > —2 or so. 

indicates that the improvement of F^^{L = U^) to the variational result is substantial when 
—3.02 < g < 0. Furthermore, the approximation of Fg up to the third order mimics the exact 
Fg better than F^\L = L^^), and substantially improves the variational result for the range 
of —6 < q < —3.02 where F^\L = L^^) is invalid. For the cases of both q < —6 and q > 0, 
F^\L = L^^) and/or F^'\L = L^^^) only slightly improve the variational result. (By the 
way, because Fg and its various approximations here vary slowly with q near their maxima, 
the humps of those curves in Fig. 6 are not evident.) 

For further and more clear illustration and comparison, we also consider the dependence of 
the various approximations here upon t. In Figs. 7, 8, 9 and 10, Fg, F^{L = L^), F^^{L = L^^) 
and F^^^{L = L^^^) are depicted for some typical values of q as the solid, short-dashed, dot- 
dashed and dotted curves, respectively, and, for the sake of clearness, we redraw them for 
the range of < t < 2 in the image. 

Fig. 7 is drawn for q = —7, and so there is not a dot-dashed curve owing to the invalidness 
of F^\L = L^^) for this case. In this figure, the long-dashed curve is F^^^{L = L^) which 
goes up when t > 5 or so, and is not given in Fig.8, 9 and 10 because it behaves too badly. 
Fig. 7 indicates that although F^\L = L^^) does not exist, F^^\L = U^^) substantially 
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FIG. 7: The dependence of Fg and its approximations up to various orders upon t at q = —7 . The 
soHd, short-dashed and dotted curves are the exact Fg and its approximations Fg{L = L ) and 
Fg {L = L ), respectively. The long-dashed curve is Fg {L = L ), going up when t > 5 or so. 
They are redrawn for the range of < t < 2 in the image. 
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FIG. 8: Similar to Fig.7, but q = -2.5, and the dot-dashed is F^^{L = L^^). 

improves the variational result and provide a better approximation for Fg. Fig. 8 is depicted 
for q = —2.5. In this figure, although the dot-dashed and dotted curves are lower than the 
solid curve, they indicate that in this case, the approximations up to the first, second and 
third orders are distinct from each other, and evidently the second- and third-order approx- 
imations more closely approach the exact value than the variational result. Fig^ and 10 are 
drawn for q = 0.5 and 1.5, respectively, which was considered in Fig.l in Ref. J5|- For these 
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FIG. 9: Similar to Fig.7, but q = 0.5 
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FIG. 10: Similar to Fig.7, but q = 1.5 

cases, the short-dashed, dot-dashed and dotted curves almost coincide and suggests that the 
approximations up to the second and third orders provide only a very small corrections to 
the variational results. 

V. CONCLUSION 

By considering the GFE for a system, this paper proposed a VPA scheme for the gener- 
alized statistical mechanics based on the Tsallis entropy. For approximating the GFE, we 
derived the truncated expressions for Fq up to the third order in the variational perturbation 
expansion, and the classical harmonic oscillator was considered in detail for an illustration. 
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The model investigation, albeit being a little academic [12| , illustrates that the approxima- 
tion up to the first order amounts to a variational method and covers the variational method 
in Ref. [5|, and the approximations up to the second and third orders improve the variational 
result and tend to approach the exact result. 

Frankly, the variational perturbation expansion technique is formally similar to the per- 
turbation expansion in Ref. jJ5|, the work in the present paper is to introduce the variational 
perturbation idea into the perturbation expansion and use the PMS for determining the 
auxiliary parameter in the VPA scheme. It is these revisions that make the variational 
perturbational approximation method be non-perturbational, take the variational result as 
the first-order approximation and systematically improve the variational result. We believe 
that the investigation in the present paper is useful, at least, for calculating the generalized 
thermodynamical functions based on the Tsallis statistics. Finally, we intend to point out 
that since path-integral formalism has been developed [1[ , it is worth while developing VPA 
method within the formalism of the generalized thermodynamics. 
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